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ABSTRACT 

The pulsar current, in the P-P plane where P is the pulsar period and P is the period 
derivative, is used as a supposedly "model free" way to estimate the pulsar birthrate 
from statistical data on pulsars. We reconsider the derivation of the kinetic equation 
on which this is based, and argue that the interpretation of the pulsar current is 
strongly model dependent, being sensitive to the form of the assumed evolution law 
for pulsars. We discuss the case where the trajectory of a pulsar is assumed to be 
of the form P = KP 2 ~ n with K and n constant, and show that (except for n = 2) 
one needs to introduce a pseudo source term in order to infer the birthrate from the 
pulsar current. We illustrate the effect of this pseudo source term using pulsar data to 
estimate the birthrate for different choices of n. We define and discuss an alternative 
"potential" class of evolution laws for which this complication is avoided due to the 
pseudo source term being identically zero. 
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1 INTRODUCTION 

The pulsar current is used to derive constraints on the birth 
of pulsars from the statistical distribution, Np(P, P), of the 
population of pulsars in the P-P plane (Figure [T]), where 
P is the period and P is the period derivative. Individ- 
ual pulsars move along trajectories on the P-P plane, be- 
coming observable ( "birth" ) at some point and disappearing 
from observation ("death") at some other point. In a steady 
state there should be a constant flow of pulsars in the P— 
P plane from the region where they are born to the region 
where they die. This flow is referred to as the pulsar current 
(iPhinnev fc Blandfordl Il98ll ; IVivekanand fc Naravanl Il98ll ; 
iDeshpande et all 19951 ). which is defined as the i ntegral over 
P of J P = PNp{P,P). In a standard model (|Lvne et all 
1 19851 ). all pulsars are thought to be born in supernovae with 
short initial periods. The expected pulsar current should in- 
crease as a function of P, reaching a plateau at P = Pq, 
with Pq corresponding to the maximum period at birth, fol- 
lowed by decline at large P as pulsars die. The birthrate is 
proportional to the height of the plateau. An early anal- 
ysis using the pulsar current technique indicated a step- 
like increase in the pulsar cu rrent at around P = 0.5 s 
|Vivekanand fc Naravanl Il98ll ). implying significant injec- 
tion at relatively large P, and also a t relatively hig h mag - 
netic fields (B > 10 12 G). However. lLorimer et all (|l993l ) 
found no strong evidence for inj ection at such relatively long 
periods. A more recent analysis (|Vranesevic et ai1l2004l ) sug- 
gests that up to 40% of all pulsars are born with periods 
in the range 0.1-0.5 s, and that over half the pulsars are 



born with magnetic fie lds > 2.5 x 10 12 G. From the latest 
pulsar current analysis, iKeane fc Kramer! (|2008l ) found that 
the inferred birthrate is not obviously consistent with the 
supernova rate. 

In this paper we reconsider the formal basis for the pul- 
sar current technique, questioning the view that it is "model 
free." The formal basis for the pulsar current technique can 
be understood in terms of an analogy between the evolution 
of pulsars and th e evolution of a one dimens ional (ID) dy- 
namical system (|Phinnev fc B landford 198 J). The period, 
P, plays the role of a generalized coordinate, and the pe- 
riod derivative, P, plays the role of a generalized velocity. 
The evolution of pulsars may be described by a force law, 
which can be written as P = f(P, P), and the solution of 
this equation gives the trajectory (the path of an individ- 
ual pulsar in the P-P plane), which depends on the specific 
function f(P, P), and on an arbitrary constant of integration 
that can be different for different pulsars . Some examples 
were given bv lPhinnev fc Blandfordl l|l98ll ). Here we choose 
a widely favoured ( "constant-n" ) model, in which pulsars are 
assumed to evolve such that they have a constant braking 
index, n, to illustrate our arguments. 

The conventional interpretation of the pulsar current 
relies on two assumptions. One assumption is that all pul- 
sars evolve according to the same law, and that this law 
depends only on P and P. The requirement that the evo- 
lution depend only on P and P implies that it must not 
depend on any additional pulsar variable, such as a chang- 
ing magnetic field, B, or a changing angle, a, between the 
rotation and magnetic axes. This requires that pulsars with 
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Figure 1. A P—P diagram showing the location of different pul- 
sars populations: normal pulsars, millisecond (together with bi- 
nary) pulsars and AXPs (anomalous X-ray pulsars) & SGRs (soft 
7-ray repeaters). Dots correspond to all radio pulsars taken from 
the ATNF pulsar catalogue. Binary pulsars are marked with cir- 
cles, AXPs & SGRs are marked with stars. Dashed arrows are 
showing direction of pulsars motion on the diagram with differ- 
ent braking index values. Solid arrows indicate the current spin 
evolution of six pulsars for which the braking index has been 
measured (see Table [TJ. Lines of constant characteristic age and 
surface magnetic dipole field strengths are shown. 

the same P and P must have the same P and the same n. 
This assumption seems inconsistent with glitching, involv- 
ing discontinuous changes in P, and other non-systematic 
changes observed in pulsar timing. The second assumption 
relates to the interpretation of P-P space as a phase space, 
which underlies the concept of a current in this space. This 
assumption is satisfied only for specific classes of force laws, 
which includes all "potential" force laws, in which the force 
f(P,P) does not depend on P. This assumption is violated 
for the constant-n model. This inconsistency can be cor- 
rected by introducing a pseudo source term, that mimics 
births or deaths, but actually corresponds to convergence 
or divergence of trajectories, due to the wrong phase space 
being used. In this sense, the interpretation of the pulsar cur- 
rent is definitely not model free: for the constant-n model 
one needs to include the pseudo source term that depends 
on n. One of our objectives in this paper is to show how this 
pseudo source term affects the interpretation of a specific 
set of pulsar data. 

Despite these seemingly serious and perhaps fatal flaws, 
it is generally accepted that the pulsar current analysis leads 
to results for pulsar birthrates that are not inconsistent with 
other data, suggesting that the pulsar current approach is 
rather robust. It seems reasonable to assume that there is an 
average universal law for pulsar evolution, on which glitch- 
ing and other random changes is superimposed. Assuming 
there is a universal evolution law, its form is not known. A 
constant-n model is one of many different models that one 
can choose, and it has the unsatisfactory feature that it re- 
quires one to introduce the pseudo source term. There are 
other choices, including "potential" models, for which this 



complication does not occur. This leads to an obviously un- 
satisfactory situation in that the interpretation of the pulsar 
current seems to depend on an almost arbitrary assumption 
concerning the mathematical form of a putative evolution 
law. After demonstrating the complications that a constant- 
n law introduces, we argue that it should be abandoned. If 
there is an (average) evolution law for pulsars, the simplest 
assumption one can make about this law is that it has a 
potential form. 

In Section[2]the kinetic equation for the pulsar distribu- 
tion is defined and discussed, and the pseudo source term is 
identified for a general class of evolution laws. In Section [3] 
we show how the inclusion of the pseudo source term affects 
results on the birthrate derived from a statistical analysis of 
pulsar data. In Section 0] we discuss existing mathematical 
models for the evolution, we extend the formal theory by 
introducing a Lagrangian formalism, and we carry out an 
analysis for a general form that includes both the constant- 
n and potential laws as special cases. We discuss the results 
in Section [5] 



2 KINETIC EQUATION FOR PULSARS 

The distribution function for pulsars (number per unit area 
of the P-P plane), Np(P, P, t), regarded as a function of the 
variables P, P, obeys a kinetic equation obtained by differ- 
entiating with respect to time. However, the resulting kinetic 
equation is not in the form (similar to Liouville's equation 
in statistical mechanics) that underlies the interpretation of 
the pulsar current. In this section we show that two forms 
of the kinetic equation differ by a pseudo source term. 



2.1 Evolution of the distribution of pulsars 

It was argued bv lPhinnev fc Blandfordl l|l98lf ) that the mo- 
tion of the observable pulsars in the P-P plane (Figure [TJ 
can be interpreted as a "fluid" flow, where the fluid den- 
sity corresponds to the density, Np(P, P,t), of pulsars in 
the P-P plane. Specifically, Np(P,P,i)dPdP pulsars with 
P between P and P+dP, and with P between P and P+dP. 
The rate of change of Np(P, P, t) with time is due to evolu- 
tion of individual pulsars and to birth and death of pulsars. 
These are described by the three terms in 

^N P (P,P) = S b (P,P)-S d (P,P), (1) 

where Sb(P,P) and Sd(P, P) are the rates, per unit time 
and per unit area of the P-P plane, that pulsars are born 
and die, respectively. 

According to the rules of partial differentiation, the time 
derivative of Np (P, P, t) on the left hand side of ([TJ) is given 
by 

p dNp(P,P,t) 
dP 

We are interested in a steady-state distribution of pulsars, 
dNp(P, P,t) /dt = 0, in which evolution and the birth and 
death rates are in balance. 
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2.2 Pulsar current 

The flow of pulsars through a vertical line in the P-P plane 
is interpreted as the pulsar current. The pulsar current is 
conventionally identified based on the kinetic equation 



d[PN P (P,P)] d[PN P (P,P)] 



dP 



+ 



dP 



S b {P,P)~S d {P,P). (3) 



On integrating ([3| over P, the integral over the term in- 
volving the P-derivate is assumed to integrate to zero. The 
integral over P of 



Jp(P,P) = PN P (P,P), 



(4) 



is interpreted as the pulsar current. (Jp(P,P) is often re- 
ferred to as the pulsar current, but it is actually the current 
density in a localized range of P.) The integral over P of 



S b (P,P)-S d (P,P) 



(5) 



dJ P (P,P) 
dP 

then expresses conservation of pulsars. This allows one to 
infer properties of the source term, Sb(P, P) describing birth 
of pulsars, and the sink term, Sd(P,P) describing death of 
pulsars, from the statistical distribution of pulsars on the 
P-P plane. 
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Figure 2. Updated standard pulsar current using 1438 non- 
recycled puls ars from Parkes surve ys. The NE2001 Galactic elec- 
tron d ensity llCordes &; La2idl2002h and the lTauris fc Manchester] 
beaming models were adopted, for sources above a 1.4- 
GMz radio threshold of O.lmjykpc 2 . Under steady-state condi- 
tions, the maximum value of current is the birthrate of pulsars. 



2.3 Pseudo source term 

The kinetic equation ([3]) differs from that obtained by com- 
bining ((2j and |T} in that the factors P and P are inside 
the respective derivatives. To derive ((3| from ([2]) one needs 
to move P and P inside the respective derivatives. Assum- 
ing that P and P are the independent variables, which is 
implicit in the use of the P-P plane, the first of these 
is trivial. For the evo l ution law P = f(P, P) assumed by 
IPhinnev fc Blandfordl (|l98ll ). the second introduces an ex- 
tra term, which we move to the right hand side and interpret 
as a pseudo source term, 



S P s(P,P) 



dfjP, P) 
dP 



Np(P,P). 



(6) 



The right hand side of ([5| is then replaced by Sb(P,P) — 
Sd(P,P) — S P s(P,P)- In order to interpret information on 
the birthrate of pulsars from statistical data, one needs to 
calculate the additional term ([6]). We discuss the effect of 
this term on the interpretation of the pulsar current for evo- 
lution at fixed n in the next section. We then argue that 
there is little evidence for the actual form of the pulsar evo- 
lution equation, that any model involves an ad hoc assump- 
tion, and that a sensible criterion to adopt in choosing an 
evolution equation is that the pseudo force term is absent. 



3 OBSERVATIONAL RESULTS 

The effect of the pseudo source term on the interpretation of 
the pulsar current is explored in this section for the widely 
favoured model of evolution at fixed n. This corresponds to 
f(P,P) = (2-n)P 2 /P in ©. The main qualitative point 
we make is that the use of the pulsar current is model- 
dependent, with the explicit dependence on the model ap- 
pearing through the pseudo-source term. 
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Figure 3. Standard pulsar current (represented with hatched fill- 
histogram-area) versus revisited pulsar currents (represented with 
the outlined and cross-hatched fill- histogram-areas). Errors have 
been excluded for clear comparison of three current histograms. 



3.1 Pseudo source term for different values of n 

We revisited our st andard pulsar current analysis 
(IVranesevic et al. I |2004 ) by including the pseudo-source 
term in our calculation. In the following analysis we 
used 1438 non-recycled pulsars that were predomi- 
nantly discovered by the Parkes pulsar surveys, of 
which all are obtained from public pulsar catalogue at 
http : //www. atnf . csiro . au/research/pulsar/psrcat , We 
accurately modelled the sensitivity threshold of the multi- 
beam survey and use this model to compute the pulsar 
birthrate. To simulate the uncertainties in the data, rather 
than choosing a fixed braking index, we chose n randomly 



4 N. Vranesevic and D. B. Melrose 



distributed about a mean, ranging from 1 to 5, and then 
calculated the pulsar current, with the same radio lumi- 
nosity threshold, Galactic electron density and beaming 
models as for the con ventional pulsar current approach 
I Vranesevic et al.|[2004 ). Figure [2] In Figure [3] we plot the 
standard pulsar current together with two revisited pul- 
sar currents, denoting their corresponding pulsar birthrates. 
There are substantial differences in currents corresponding 
to different values in n. For n > 2, the S ps {P, P) is negative, 
which reduces the pulsar current, implying that the neglect 
of the pseudo source term leads to an overestimate of the 
current. For n < 2, S ps (P,P) is positive, implying that the 
neglect of the pseudo source term leads to an underestimate 
of the current. 

Our simulations for different runs of braking index are 
shown on Figures [4] a) , b) & c). Each plot displays two his- 
tograms, the revisited pulsar current on the left, and its cor- 
responding histogram of randomly chosen braking indices 
on the right. Matching birthrate and average braking in- 
dex with its standard deviation are displayed on each his- 
togram. The effect of the assumed value of n on the esti- 
mated birthrate for given data can be seen from the plots: 
n < 2 (n > 2) the inferred birthrate is higher (lower) than 
when the pseudo source term is neglected. 



3.2 Measured braking indices 

iLvnd \2004 ) and iLivingstone. Kaspi fc Gavriij (|2005h re- 
ported braking index measurements for six young pulsars 
(Crab, B0540-69, Vela, J1119-6127, B1509-58, and J1846- 
0258). Five have an unambiguous measurements of n, via 
phase-coherent timing, and these are marked by a star in 
Table[T] We show their implied trajectories in the P-P plane 
by arrows in Figure [T] 

There are only a few pulsars that are potential candi- 
dates for accurate measurement of n, which requires that 
they satisfy the following. Firstly they must spin down 
sufficiently quickly to allow a useful measurement of P. 
Secondly, the position of the pulsar should be accurately 
known at the ~1 arcsecond level. Thirdly, the spin-down 
must not be seriously affected by glitches, sudden spin-ups 
of the pulsar, or timing noise. Typically, glitches begin to 
seriously affe ct smooth spin down at characteristic ages o f 
~ 5-10kyr (|McKenna fe Lvne!ll990l ; lMarshall et alj|2004h . 
Thus many of the pulsars that may spin down fast enough 
for a me asurement of n are ir r etrievably contam inated by 
glitches (|Shemar fc Lvne!ll996l ; IWang et al.ll200ll ). Timing 
noise varies fro m object to object, is rough ly correlated with 
spin-down rate (|Arzoumanian et all 1 19941 ) and can prevent 
a measurement of n in a finite data set in an unpredictable 
way. All six trajectories have n < 3, three of them have up- 
ward slopes (n < 2), and none of them is directed towards 
the main body of pulsars. These data do not support the 
constant-n model, which requires the same value of n for all 
pulsars at all P, P. 



4 ALTERNATIVE EVOLUTION LAWS 

For the pulsar current to be a well-defined concept, there 
must be a universal law of the form P = f(P, P) obeyed 
by all pulsars. This is clearly not the case. For the pulsar 
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Figure 4. Revisited pulsar current for different random braking 
indices. Each of a), b), and c) plots represents combination of pul- 
sar current histogram on left (with birthrate written at the right 
top corner of the plot) and its corresponding histogram of ran- 
domly chosen braking indices on right (with the average braking 
index and its standard deviation written at the right top corner 
of the plot). 
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Table 1. Braking index measurements for six pulsars. Also given are the pulsar period, period derivative, period 
second derivative, and characteristic age. 

Pulsar names P P P t c braking index 

J2000; B1950 (ms) (10~ 15 ) (s^ 1 ) ( yr) n 

Crab* 

J0534+2200; B0531+21 33.085 423 -3.61 x lCT 24 1,240 2.51(1)^ 

J0540-6919; B0540-69* 50.499 479 -1.6 x 10~ 24 1,670 2.140(9^, 
Vela 

J0835-4510; B0833-45 89.328 125 11,300 l-4(2£ 

J1119-6127* 407.746 4022 -8.8 x 10~ 24 1,160 2.91(5}£ 

J1513-5908; B1509-58* 150.658 1540 -1.312 x 10~ 23 1,550 2.839(3)^ 

J1846-0258* 325.684 7083 728 2.65(1]/, 



Demiahski & Proszviiskil jl983l);lLvne et al.l l| 19881 . 1 1993h 



Livingstone. Kaspi fe Gavrii I J2005h 
Lvne et al.l l|l996MDodson et al.1 l|2002h 
Camilo et all ll2000h 

Manchester et al.l dl985h ; iKaspi et alj ll994) 



Livingstone et al.l l l2006h 



current to be a useful concept, such an evolution law must 
apply at least in some average sense for all pulsars. The 
constant-n evolution law assumed in the previous section, 
although widely favoured, is essentially ad hoc. In this sec- 
tion we introduce a general (Lagrangian) formalism, that 
assumes that there is a unique evolution law but makes no 
assumptions on its form. We use this formalism to identify 
a simpler class of evolution models that we refer to as po- 
tential models. 



radiation has a frequency equal to the rotation frequency 
of the pulsar, well below the plasma frequency in the pul- 
sar wind and in the ISM, and it cannot escape to infinity, 
as the model requires. The model predicts that the mag- 
netic and rotation axes align on the spin-down time scale. 
There is a theoretical argum ent that the str ucture of the 
star may slow the alignment (| Goldreichl 1 1 970h . and there is 
evidence for older pulsars hav ing smaller of values of sin a 
l|Weltevrede fc Johns ton 2008), but the prediction of the 
model itself is inconsistent with pulsar data. 



4.1 Specific models 

Before introducing a general theory for a possible pulsar 
evolution model, we summarize the models that underlie 
the discussion in the previous section. 



4-1.1 Vacuum dipole model 

The vacuum dipole model has a contradictory status in pul- 
sar physics. One the one hand, it is the basis of our inter- 
pretation of pulsars, notably Bsina oc (PP) 1 ^ 2 , and it is 
the only model for which the evolution is well determined. 
On the other hand, the model is clearly inapplicable to ac- 
tual pulsars, and the actual evolution is inconsistent with 
observations. The constant n — 3 model is often assumed to 
describe evolution in the vacuum dipole, but this is incor- 
rect. 

The model is that a (structureless, point) rotating mag- 
netic dipole in vacuo radiates electromagnetic radiation that 
escapes to infinity, carrying off energy and angular momen- 
tum. The loss of energy causes P to increase on a character- 
istic spin-down timescale, and the torque causes sin a to de- 
crease, implying that the rotation and ma g netic axes tend to 
align on the spin-down time scale iMichell (11991! ). The brak- 



ing index is n — 3 + 2 cot a (Michel 1991), which is always 



greater than 3, and which increases as sin a decreases. 

This vacuum-dipole model, at least in its simplest form, 
is unacceptable for pulsars. The putative electromagnetic 



4.1.2 Then = 3 model 

The constant n = 3 model is often confused with the vac- 
uum dipole model. In the n = 3 model, it is assumed that the 
evolution of pulsars is at constant Bsina, and this is com- 
bined with the vacuum dipole model so that the assumption 
becomes evolution at constant PP. Note that this is incon- 
sistent with the actual vacuum dipole model, which implies 
systematically decreasing sin a. Thus, the basic assumption 
in the constant n = 3 model is that the trajectory of pulsars 
in the P-P is of the form P — K/P, where K is a constant. 
The force law implied by this relation is P = —P 2 /P and 
this implies a braking index 



n = 2 - PP/P 2 



(7) 



equal to 3. Although a vacuum model modified by assuming 
Bsina remains constant is widely favoured it is not sup- 
ported by any detailed theory for the evolution. 



4-1.3 n = constant model 

Another model that is used widely is based on the assump- 
tion that the braking index is constant but not necessarily 
equal to 3; we refer to this as the n = constant model. Ac- 
cording to (0, n = constant implies P = (2 — n)P 2 /P, 
and integrating this equation leads to the trajectory P — 
KP 2 ~ n , where K is a constant. This is the model used for 
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illustrative purposes in the previous section. However, the 
assumption n = constant on which it is based is also not sup- 
ported by any detailed theory for the evolution, and should 
be regarded as ad hoc. 



4-1-4 Potential models 

Granted that one needs to make an ad hoc assumption in 
choosing a possible universal law governing pulsar evolution 
in P-P, it is desirable to choose a law that does not intro- 
duce unnecessary complications. As shown in detail above, 
the n = constant model leads to the need to include a pseudo 
source term in the interpretation of the pulsar current. A 
model that avoids this complication is one in which P de- 
pends only on P. We refer to this as a potential model. The 
law is written in the form P = —C'(P), where C(P) is the 
potential in the theory developed below. This theory is also 
not supported by any detailed theory for the evolution, and 
is also ad hoc. However, a potential model has the advantage 
of being the simplest model, and the best justification for it 
is that one should explored the simplest model first. 



4.2 Lagrangian formulation 

In a Lagrangian formulation of dynamics the independent 
variables are the generalized coordinates and generalized 
velocities. In the present case, the generalized coordinate 
is identified as P and the generalized velocity as P. A La- 
grangian, L(P,P) may be identified by requiring that the 
Eulcr-Lagrange equation 



d dL(P, P) dL(P, P) 



dt dP 



dP 



0. 



(8) 



reproduce the equation of motion. 

Consider the evolution law ((Phinnev fc Blandfordfl98lf ) 



(2 - n)P 2 



P 

T 



(9) 



This has a Lagrangian counterpart for r — > oo. The relevant 
Lagrangian is a particular case of 

L(P,P) = \A(P)P 2 + B(P)P-C(P), (10) 
for which Q implies 

2A(P) A(P)' ( ' 

where the prime denotes a derivative in P. The form (O is 
reproduced for A(P) = p- 2 ( 2 "™), C(P) = 0. 

The trajectory of a pulsar in P-P space is found by 
solving the equation of motion. This is facilitated by noting 
the analogy with a one-dimensional particle in a potential, 
with A(P) playing the role of the mass and C(P) playing 
the role of the potential energy. The energy-integral (sum of 
kinetic and potential energies equals a constant) gives 



\A{P)P 2 + C(P) = E, 



(12) 



where E is a constant of integration (the counterpart of the 
energy). A second integral gives 



dP 



' A{P) 
E-C(P) 



1/2 



2 1/2 t. 



(13) 



The example considered bv lPhinnev fc Blandfordl (| 198 if ) is 
reproduced by writing E = \v 2 , x = vt. 



4.3 Potential models 

There is no pseudo source term in the kinetic equation for a 
Lagrangian of the form (j 10|) with A(P) = 1. The generalized 
momentum, dL/dP, is then equal to the generalized veloc- 
ity, so that the P-P is a phase space in the sense required 
for Liouville's theorem to apply. In this case pip gives 



P = -C'{P), 



(14) 



and 1(13) 1 apply with A(P) = 1. The braking index in 

this model is 



PP 
P 2 " 



PC'(P) 
2{E-C(P)] 



(15) 



A simple limiting case is E 3> C(P), when the evolution law 
reduces to P ~ constant, with a braking index n ~ 2. 

In these two examples of evolution laws, the trajectories 
of individual pulsars depend on a constant of integration, 
which depends on the particular pulsar. For trajectories of 
the form P = KP 2 ~" this constant is K, and for the po- 
tential law 1(120 . the constant is E. For any given law, the 
constant of integration is determined for a particular pulsar 
by the values of P and P at birth, and the subsequent evo- 
lution is determined by (|12|) with E fixed, but different for 
different pulsars. 



5 DISCUSSION AND CONCLUSIONS 

The theory for the pulsar current is based on the assumption 
that there is a universal evolution law for pulsars, and that it 
depends only on P and P. This requirement is not satisfied, 
with glitching and other random changes implying that the 
trajectory of at least some pulsars on the P-P plane is in- 
fluenced by other effects. Nevertheless, it is possible that the 
evolution of all pulsars does follow a universal law in some 
average sense, and this would suffice to justify the pulsar 
current technique. Assuming there is an average universal 
evolution law that applies to all pulsar, it is clear that we 
have little information on what this law is. A constant-n law 
has been used widely in the literature. We assume this law 
in Section[3J argue that such a law invalidates a conventional 
pulsar current analysis, and show how the error can be cor- 
rected by introducing a pseudo source term. We argue that 
there is no theoretical basis for a constant-n law, and that 
the complication it introduces through pseudo source term 
makes it an undesirable choice for the form of the evolution 
law. We argue that the choice of a potential law for the evo- 
lution avoids the complication of a model-dependent pseudo 
source term. 

We discuss evolution laws in Section 2) We comment on 
the misconception that a constant n = 3 law corresponds to 
the vacuum dipole model, whereas it actually corresponds 
to a hybrid model, with Bsina oc (PP) 1 ' 2 determined by 
the energy loss rate in the vacuum dipole model, and with 
Bsina assumed constant. The actual vacuum dipole model 
implies that sin a decreases (alignment of the magnetic and 
rotational axes) on the spin-down time scale, and this does 
not occur. The vacuum dipole model is not applicable to 
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actual pulsars, for this and other reasons. A general class of 
evolution models is developed in Section[4] this class includes 
the constant n model, with P — (2 — n)P 2 /P, and potential 
models, with P = ~C'(P) where C(P) is the potential. The 
evolution law is integrated to find the trajectory of the pulsar 
in P—P space. The trajectory can be written P = KP 2 ~ n 
for a constant-n model, and as P = \/2[E — C(P)] 1 / 2 for a 
potential model. The trajectory depends on a constant of in- 
tegration, which is different for different pulsars. A physical 
model for the evolution is required to give a physical inter- 
pretation to such constants of integration. In the constant 
n = 3 model, the constant K is proportional to Psina, 
but there is no obvious interpretation of K for other models 
with constant n / 3. A physical basis for a specific potential 
model is required to give a physical interpretation to E. 

Assuming a potential law opens up a new way of think- 
ing about the evolution of pulsars in the P—P plane. In our 
formulation the potential, C(P), is an arbitrary function of 
P. By exploring different choices of C(P) one can ask what 
choice (if any) best fits the data. Although we do not at- 
tempt to do this here, we speculate briefly on one implica- 
tion. Suppose that the "potential energy" is much smaller 
than the "kinetic energy" [C(P)| -C \P 2 ■ Then all trajecto- 
ries have P ~ constant, implying a braking index n ~ 2. The 
measured values of n include n < 2 and n > 2, which could 
be interpreted either according to (|15p . requiring C'(P) < 
and C'(P) > 0, respectively, or interpreted as random vari- 
ations about an average universal n ~ 2 for young pulsars. 
Perhaps average braking indices can be estimated for a large 
number of pulsars over a wi de range of P, as suggested by 
Ijohnston fc Gallowavl (jl999h for example. If such an estima- 
tion of n shows some systematic trend with P, (|15p would 
provide some constraint on the possible form for C(P). With 
such data one should be able to draw conclusions on whether 
simple examples potential law are consistent with observa- 
tions; simple models include C(P) oc P a with a > 0, and 
C(P) oc -P~ h with b > 0. If there is no systematic trend, 
even in an average sense, this would cast serious doubt on 
the validity of the pulsar current technique, and conclusions 
drawn using it. 

In conclusion, the pulsar current approach is not "model 
free" and its interpretation depends on assumptions made 
about the form of the evolution law for pulsars. The theory 
depends on this law being universal, and although this is 
clearly not the case in an exact sense, it may be the case in 
some average sense. In modelling such an average law, the 
constant-n assumption should be abandoned, and a "poten- 
tial" law should be the preferred choice on which to base 
future discussions. 
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